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QCD sum rules are evaluated at finite nucleon densities and temperatures to determine the change 
of mass parameters for the lightest vector mesons p, uj and <fr in a strongly interacting medium. For 
conditions relevant for the starting experiments at HADES we find that the in-medium mass shifts 
of the p and u> mesons are governed, within the Borel QCD sum rule approach, by the density and 
temperature dependence of the four-quark condensate. In particular, the variation of the strength 
of the density dependence of the four-quark condensate reflects directly the decreasing mass of the 
p meson and can lead to a change of the sign of the ui meson mass shift as a function of the density. 
In contrast, the in-medium mass of the <f> meson is directly related to the chiral strange quark 
condensate which seems correspondingly accessible. 



The in- medium modifications of the light vector mesons (p, w, </>) receive growing attention both from theoretical 
and experimental sides. On the theoretical side there are various indications concerning an important sensitivity 
of vector mesons to partial restoration of chiral symmetry in a hot and dense nuclear medium. In particular, at 
finite temperature the vector and axial-vector correlation functions, which are related to the meson spectral densities, 
become mixed in accordance with in-medium Weinberg sum rules |l| . At low temperature this mixing can be expressed 
directly via vacuum correlators in a model independent way Additionally, as shown within lattice QCD || and 
various effective model calculations [Q, the chiral quark condensate as order parameter decreases with increasing 
temperature and baryon density. Within the QCD sum rule approach, considerable in-medium modifications of 
vector meson masses even at normal nuclear matter density have been predicted m . 

On the experimental side there is the idea to probe the in-medium modifications of vector mesons, in particular 
mass shifts, by measuring dileptons (e + e~) from meson decays. This is the primary motivation of the presently 
starting experiments with the HADES detector at the heavy-ion synchrotron SIS at GSI Darmstadt. At higher 
beam energies the heavy-ion experiments of the CERES collaboration |^| evidenced already hints to a noticeable 
modification of the dilepton spectrum which can be reproduced under the assumption of a strong melting of the p 
meson in a dense, strongly interacting medium at temperatures close to the chiral transition P,^0[. 

While numerous evaluations of QCD sum rules have been performed during the last decade, the majority of them 
deal either with cold nucleon matter or a hot pion medium. At the same time, to extract the wanted information 
on the behavior of in-medium mesons via measurements in heavy-ion collisions at SIS energies one has to study the 
case of finite baryon density and temperature. In this paper we present evaluations of the QCD sum rules for the 
light in-medium vector mesons in a density - temperature region which is relevant for experiments like the ones at 
HADES. We study systematically here the relative numerical contributions of the four-quark condensate to the QCD 
Borel sum rule to find out the general trends of the vector meson mass dependence on density and temperature with 
respect to variations of the poorly known four-quark condensate. 

Our paper is organized as follows. In section 2 we recapitulate the general structure of the QCD Borel sum rule. 
The operator product expansion and evaluation of the needed condensates are summarized in section 3. In section 4 
we discuss the form of the hadronic spectral density. The numerical evaluation of the sum rule is presented in section 
5. The summary can be found in section 6. 
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I. INTRODUCTION 



II. QCD SUM RULES AT FINITE DENSITY AND TEMPERATURE 



Following the approach developed in |Tl],|T2| we employ for the QCD sum rules (QSR) at finite nucleon chemical 
potential piN and temperature T the retarded current-current correlation function in a medium 




(1) 
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where x = (x°, x), q = {q°,q) and TZJ l _ l (x)J u (0) = Q(x°)[J^(x) 7 J v (0)] with the conserved vector current J M ; (■ ■ -) mw ,t 
denotes the grand canonical ensemble average. 

We consider vector currents of QCD with isospin quantum numbers of the respective vector mesons specified 
as = Mu'YuU =F drffid), where the negative (positive) sign is for the p (u>) meson and = s7 M s is for the <j) 
meson current. Since we focus on the limit q —> in the rest frame of the medium, only the longitudinal invariant 
Ilf = n- RM Al /(— 3q 2 )\^o of the correlator ([!]) is needed in our analysis. 

Due to the analyticity in the upper half of the complex q° plane the retarded correlation function Hf; satisfies the 
standard dispersion relation in a medium 

n^^^lf^J*^, (2) 

Ti" Jo s-{q + ie) 2 

where a subtraction has been omitted. For large Q 2 = —q$ > one can evaluate the correlator II^ by expanding the 
product of currents in (Q) by means of the operator product expansion (OPE). The result can be written as 

oo „ 

n£(Q 2 ) = -C lnQ 2 + ]T^, (3) 

n— 1 

where the quantities C n include the usual Wilson coefficients and the condensates as well. It is remarkable that in 
the considered region (q° ^ real) the OPE for Il£ is the same as the OPE for the causal (Feynman) correlator. This 
can be used to get an explicit expression of (|3j) in a simple manner. 

As usual within the QSR the spectral density phad(s; pn,T) = ilmll^ (s; pjy, T) is modeled by several phenomeno- 
logical parameters related in particular to the forward scattering amplitudes of the external current and the 
constituents of the medium. 

Performing the Borel transformation with appropriate mass parameter M 2 fl3|| of the dispersion relation eq. @j 
and taking into account the OPE ([}]) one gets the basic equation for the QSR evaluation in a medium as 

(°° C \ 



o 



which formally looks similar to the corresponding equation derived in |e| for the Borel sum rules at finite temperature. 



III. OPE AND EVALUATION OF THE CONDENSATES AT FINITE p, N AND T 

The coefficients C n , which define the r.h.s. of the basic equation (Q), include the Wilson coefficients and the grand 
canonical ensemble average of the corresponding products of quark and gluon field operators. Their general structure 
up to n = 3 is given, for instance, in G3j. For p and ui mesons one has 



3m 2 



C 2 = m q (uu)^ NiT + 7^( — G )hn,t - ^-QriSTiwr^Duu))^^, (7) 



q=u,d,s 

,\&qi 1 q v q x q a 



Q 4 



{ST{u ltl D v D x D a u))^ N , T , (8) 



where again the minus (plus) sign corresponds to p (lu) meson; a s (p 2 ) is the strong coupling constant which is 
taken at the renormalization point p? = 1 GeV 2 according to common practice; m q = \{m u + md) stands for the 
average light quark mass; G 2 — G^G^ with the gluon field strength tensor G£ v . The covariant derivative is defined 
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as L> M = dp + yA*\ a with SU(3) color matrices A a normalized as Tr(A a A fc ) = 2S ab . The operator ST creates a 
symmetric and traceless expression with respect to the Lorentz indices. 

To get the /i/v and T dependence of the condensates, entering the coefficients 62,3, we assume that the system, at 
small density and temperature, can be described by non-interacting nucleons and pions. For small enough /i^r and T 
the particle gas is dilute and the grand canonical ensemble average of an operator O can be approximated by 

(0) m , T = (0)0 + 3 J {2 ^ P 2Ep n B{<p)\OU{ P )) 

/d 3 k 

where {O)o is the vacuum expectation value, and Hb = [e Ep / T — l]^ 1 and rip = \e^- Ek ~ tJ - N ^ T + are thermal Boson 
and Fermion distributions; we use the covariant normalization of the one-particle pion (nucleon) state: (tt(jj)\tt(p')) — 

(2tt) 3 2E p 6(p-p') {(N(k)\N(k')) = (2tt) 3 2E k 6 3 (k - k')), where E p = yjp 2 + ml (E k = \Jk 2 + M 2 N ) with 
(jun) as pion (nucleon) mass; the vacuum is normalized as (0 1 0) = 1. 

Let us consider first the matrix elements of the scalar (Lorentz invariant) operators. For the vacuum quark and 
gluon condensates, (qq)o with q = u,d,s and (^G 2 )o, we employ values presented in table 1. For the vacuum four- 
quark condensate we adopt the vacuum saturation hypothesis |13j] which is valid in the large- N c limit (N c is number 
of colors) so that the four-quark condensates become proportional to the squares of the quark condensates. To control 
a deviation from the vacuum saturation we introduce the parameter kq in the following manner 

16 

((wA a q) 2 }o = -((qi^ a q) 2 )o = -Y K o{qq)l (10) 

The case kq — 1 corresponds obviously to the exact vacuum saturation as used, for instance, in Jl2] |. 

The pion and nucleon matrix elements of scalar operators which contribute to (0) do not depend on the particle 
momenta and, therefore, can be calculated in the limit of vanishing momenta. The pion matrix elements of quark 
operators can be expressed via vacuum condensates by means of the soft pion theorem Jl2]| 

(7r|uu|7r) = (7r\dd\ir) = — —(uu)o, (7r|ss|7r) = (11) 

f4 

with the vacuum pion decay constant f„ — 93 MeV. 

Among four-quark operators in eq. (|^), only the term (u7 A1 75A a u — dj^5\ a d) 2 = 0$ , which contributes to the p 
meson QSR, has a non- vanishing pion matrix element. Applying the soft pion theorem twice leads to 

(n\0^\7r) = -^L{(u 1 , l5 \ a u)\. (12) 

"J J 7T 

Using the vacuum saturation assumption with the corresponding parameter kq in eq. ([lO]) one gets 

(n\O^\n) = -—^K (uu)l (13) 

The pion expectation value for the scalar gluon condensate is obtained as usual by employing the QCD trace 
anomaly of the energy-momentum tensor [n2|n3] and the Gell-Mann - Oakes - Renner relation 

H-G 2 k) 8 m 2 . (14) 

7T 9 

By a comparison of the quark condensate in a nucleon medium obtained by means of the Feynman - Hcllmann 
theorem, applied to the ground state of nuclear matter in Fermi gas approximation |l4| , 
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<«>wr-<«>o = 2d^. e = A j^ <k J^V k2 + M *> k F=\J»N- M N (15) 
and eq. @ taken at T = one can get the nucleon matrix element of the scalar quark operator as 

(NMN) = **™, q = u,d, (16) 
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where o~n 



, <iM/v / duiq is the nucleon sigma term. Similar steps for the s quark condensate result in 



(N\ss\N) = y 



M N (JN 



(17) 



with the parameter y defined via ycrjv = 2m s dMN / dm s . Our choice of the numerical values of parameters entering 
eqs. (|l6| |l7|) is given in table 1. 

To get a proper ansatz for the nucleon matrix element of the scalar four-quark condensates we extend the widely 
used ground state saturation assumption for the in-medium four-quark condensate at T = in the following way 
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((nu 7 5 A a <zr) AIN = -((qj^Xq) 2 )^ N = — K (n N )(qq) 



(18) 



where the density dependent factor k(tijv) is introduced to control a deviation from the exact ground state saturation 
in analogy to the vacuum saturation in eq. (|To|) . As pointed out in |15| the factor k(jin) reflects the contribution 
from the scalar low-energy excitations of the ground state and seems to be weekly dependent on the nucleon density 
npf, so that one can use k(jin) = const as first approximation. In this case the density dependence of the four-quark 
condensate appears only via the density dependence of the quark condensate squared. In linear density approximation 
it is given by 



= (qqfo + (qq)o(N\qq\N) 



n N 



M 



N 



(19) 



We go beyond the above first approximation and assume here also a possible linear density dependence of k{un), 
keeping in mind that still /t(njv) = kq at tin — in accordance with the vacuum case in eq. (|10|). To leading order in 
the density this assumption can be expressed as correction to the second term of the r.h.s. of eq. ( |l9| ) by a constant 
factor kn- As a result, our parameterization of the four-quark condensate at T = has the form 



<(97/x75A a g) 2 ) PiV = 



16 



o K o 



1 



k n (N\qq\N) n N 
k (qq) M 



N 



(20) 



The limit kn = kq = 2.36 is used in fL5| , while the parameterization kn — 1-4 and kq — 3.3 with (qq)o = (— 230MeV) 3 
corresponds to the choice in ||. Below we vary the parameter kn to estimate the contributions of the four-quark 
condensates to the QSR with respect to the main trends of the in-medium vector meson mass modification. 

The needed ansatz for the nucleon matrix element of the scalar four-quark condensate can be extracted then from 
the direct comparison of our parameterization in eq. ( p0| ) and the general expression (^|) for the condensates via the 
matrix elements (the latter ones to be taken at T = 0) as 
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(N\(<11^5X a q) 2 \N) = —{qq) Q (N\qq\N) k n . 



(21) 



Since the matrix element fclj) does not depend on particle momenta and temperature it can be employed also with 
the approximation in eq. (£3|) for evaluating the four-quark condensates in the general case with T ^ and [In ^= 0. 

The gluon condensate in the nucleon is obtained in the same manner as for pions by using the trace anomaly and 
eqs. (0 0) 



(iV|^G 2 |A0 

7T 



~M N M%, 



(22) 



where is the nucleon mass in the chiral limit; for numerical values see again table 1. 

The pion matrix elements of twist-2 operators, which are symmetric and traceless with respect to Lorentz indices, 
can generally be written as 



(7r(p)|STg7^A,g|7r(p)} = -i [p^Pu ~ -.Q^P A% (n ) 



for a = 



u, d, s and for mass dimension-4 operators, and 



1 



(23) 



(7r(p)|§T^A,I>A.D ff g|7r(p)) 



PtiPvPXPa 



4cS 



"g" (PfiPvgXa + P^Pxgua + P^PagXv + PuPxg^a + PuP^g^X + PxPag^u) 



AIM 2 ) 



(24) 
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for mass dimension-6 operators, respectively. The coefficients A\ v and A\ w are defined by 

A\ ^p 2 ) = 2 / dxx^ 1 [<^(x, M 2 )+^(x, M 2 )] , * = 2,4 (25) 
Jo 

where q 7r (x,p) is the quark distribution function inside the pion at the scale /x 2 . 

The nucleon matrix elements for twist-2 operators of mass dimension-4, (N(k)\STq r y IJi D l/ q\N(k)), and mass 
dimension-6, (N(k)\STq~'y l j,D l ,D\D cr q\N(k)), have the same structure as in eqs. (^3| - |25| ), but with the quark dis- 
tribution qisf(x,fj,) inside the nucleon at scale /i 2 and AJ N (p 2 ) = 2 dx a; 2-1 [qN(x,p 2 ) + Qn(x,p 2 )] . Our choice of 
the quantities A?^^ N ^ is listed in table 1. 

Using the obtained pion and nucleon matrix elements of the corresponding operators and performing the grand 
ensemble average eq. (^|) one finally gets the values of the coefficients C2,3 entering the basic equation (0) for p and 
uj mesons as 



c 2 


= 92 + .92 +0,2, 






(26) 


c 3 


— 94 + 04, 






(27) 


92 


= m q (uu) + 2M N a N YI? + -m^Tf 




(28) 


.92 


= — (^G 2 ) - — M N M a N I? - —mill 
2V 7T 70 27 w 1 18 ^ 1 




(29) 


94 




' k n AM N a N N 36 
«o m g (-(uu)o) 1 7/ 2 


fix 


(30) 


(l 2 


4 S 

— M 2 A u+d T N J- _|_ m 2 4"+<2 j-t 1 4«+<i 

— JVtfjy A 2 JV J x + g^2, TV J 2 +J m ir A 2 : TT 1 1 + A 2,TT 


is 


(31) 


O4 


- 5 M 4 A u+d T N 

— 2% 4,iV J l 


20 M 2 tN 16 rJV 








m 4 4 u+d r 71 ' Km 2 A u+d T n AA u+d T v 




(32) 



where Y = 1, = 1 for p and u; mesons, £ p = 1 for the p mesons, and elsewhere = 0. The sequence of 
replacements m q m s , (uu)o — > (ss)o and Y - — > ym s /m q , A^^ N ^ — ■> A?^^ holds for the meson. The above 
integrals are 

The density and temperature dependence of the scalar quark and gluon condensates q2 and (72 are exhibited in fig. 1 
in |jT7f . One observes a striking linear density dependence, while for small temperatures there is only a tiny change 
of the condensate. A suitable approximation for (qq) nN ,t / '(99)0, (ss} nN: T/(ss)o and (^tG 2 )„ N; t/(yG 2 )o at small 
temperatures is 1 — an/no with a w 0.32, 0.09 and 0.08, respectively (we use no = 0.17 fm~ 3 ). 



IV. HADRONIC SPECTRAL FUNCTION 



To get insight for modeling the hadronic spectral density /?had(s; /ijv, T) — ^IuiII^ (s; /ijv, T), which enters the l.h.s. 
of the basic equation (Q), one can decompose the in- medium correlator (q; /xjv,T) in the same manner as in eq. (|^) 
within the dilute gas approximation by restricting on not too large values of /xjv and T 



II* ( g ; Mjv ,T)=n« (q;0,0) 



+3 J WF^ nBT;M ® +A l lW^ nFT "» {q;k) ' (34) 

where n*„ (q; 0) is the vacuum expectation and 

^(9;p) = ijd 4 x e^(n(p)\TZJ,(x)J,(0)\TT(p)), (35) 

^(9; *) =* [d 4 x e^(N(k)\TUp(x)J u {0)\N(k)) (36) 
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are the forward scattering amplitudes of the external current J M for pion or nucleon, respectively. 

For the vacuum QSR, only the first term in eq. (|34| ) is obviously non- vanishing. The corresponding spectral density 
for a given vector meson channel was successfully modeled in [T^ ] by means of a resonance peak and a continuum 
(perturbative QCD) contribution. In a medium one has to take into account also the contributions arising from the 
imaginary parts of the second and third scattering terms in eq. j3^). To model such contributions we consider only 
absorption like processes related to the matrix elements (ttIJ^tt) and (N\ J M |iV). This is still in accordance with the 
dilute gas approximation when only the one-particle states are taken into account. In the limit q —* the above 
scattering terms, called usually Landau damping terms, can be calculated exactly (cf. ||jll]] for details). As a result, 
our parameterization of the hadronic spectral density in the vector meson channel V is given by 



v 

Phad 



V 
Phaa 



V 
Ptcs 



V 
Peon 



(37) 



with the usual splitting into the Landau term PLan 



pY es and the continuum contribution p^ on 



c e(s 



(Pscatt + Psca±t)3( s ) (specified further below), the resonance part 
s^). The continuum threshold Sq separates the resonance part 



of the spectrum from the continuum part. 

Without specifying the resonance part one can define the resonance mass parameter my as the normalized first 
moment via 



m v 



j;°d sp Y cs ( s y- s/M2 



(38) 



which is evident in the zero- width approximation where pY es = Fy S(s — m v ). One should underline that, according 
to eqs. (|| [38|), only the weighted integrals over pY es (s) enter the QSR. That means that various shapes, such like the 
pole approximation or particular Breit-Wigner parameterizations or double-peak distributions etc., can deliver the 
same value of the mass parameter my supposed the corresponding integrals over pY cs are fixed. 

The above functions p^It^ entering the Landau damping term in eq. (E?7]) of the spectral density are given for the 
p meson by 



p,N 
Pscatt 



p,7r 
Pscatt 



1 



48tt 2 
1 

24tt 2 



dui hp\ 1 — 



4Mt 



dio ns\/l 



4m? 



AM 2 N 



4m 



(39) 
(40) 



and n F = [ e ("- 2 ^«)/ 2T + l]" 1 , n B = [e w / 2T - For the u> meson 
due to the different isospin structure of the p and ui mesons |T~ 
and nucleon gas is negligible 



PscZtt = due t0 symmetry, and p"£ tt = 9p^ tt 
For the <fi meson the Landau damping for the pion 



V. EVALUATION OF QCD SUM RULE 

Inserting the hadronic spectral function eq. (|37|) in the basic equation (^), taking the derivative with respect to 
M~ 2 , and using the definition of the mass parameter my in eq. ( J38| ) one finally gets the sum rule in the form 



C (l 






- C 2 M- 4 - C 3 M- 6 


C (l-e- s VM^ . 


f CiM" 2 


f C 2 M- 4 - 


f iC 3 Af-6-(p^ t+p ^ tt )M- 2 



(41) 



with the above coefficients Co... 3 from eqs. (||, [| - ^). The identification of my with the vector meson mass is 
suggestive, however, strictly speaking only valid in the limit of narrow resonances. Then the pole approximation is 
applicable and, as shown in [Q, the width of the resonance is calculable from the pole residue Fy afterwards. Having 
this in mind we focus on the in-medium change of the mass parameter my which is associated to the vector meson 
mass. 

The quantity my{M 2 ) still needs to be averaged within a certain Borel window M m { n ■ ■ ■ M max . The continuum 
threshold Sq is determined by requiring maximum flatness of my(M 2 ) within the Borel window. We use a sliding 
Borel window determined by the "10% + 50% rule" p0| , pl| , i.e., the minimum Borel mass M^ in is fixed such that the 
terms of order 0(1/M 6 ) on the OPE side contribute not more than 10%, while the maximum Borel mass M^ ax is 



G 



determined by the requirement that the continuum part equals to the resonance contribution in the spectral function. 
(A fixed Borel window with M m i n = 0.8 GeV 2 and M max = 1.5 GeV 2 (2.5 GeV 2 for the <j> meson) delivers the same 
results.) 

The results of our QSR evaluation for the density and temperature dependence of the resonance mass parameter 
uiy are exhibited in figs, [l] - ||. As seen in figs. 0a and b the approximately linear dropping of m p with increasing 
density appears to be stable with respect to variations of the temperature in wide region < T < 140 MeV. In 
particular, if one parameterizes the four-quark condensate (see eq. (2^)) by kn — 1-4 and kq = 3.3 according to ||, 
the density dependence can be approximated by m p = m? p (l — a p npf /no) with a p w 0.16 remaining constant with 
respect to variations of the temperature within 30% accuracy. In fig. m a one also observes that at fixed density the 
mass parameter m p suffers only small changes when increasing the temperature up to 100 MeV. This is in line with 
the Eletsky-Ioffe mixing theorem on the in-medium behavior of the vector and axial-vector correlation functions in 
the chiral limit P, p2| , p3[| ■ Up to order T 2 , at tin = the OPE eq. (^) for the p meson correlator in a pion medium 
can be rewritten in the following form 

nf(Q 2 ;T) = [n p (Q 2 ;0) + e(n ai (Q 2 ;0)-n p (Q 2 ;0))] , (42) 

where n pai (Q 2 ;0) = IT^ (Q 2 ; 0) are the vacuum vector and axial-vector correlators and e = ^-.B (^) with 
B(x) = ^ J°° dy\Jy 2 — x 2 /(exp(y) — 1) (note that B(Q) = 1). In deriving eq. ( f42|) we use the dilute gas averaging 
@ and the vacuum saturation hypothesis ( [To| ) which obviously give IT Ol (0) — n p (0) = 32 q" s (qq)oK . From eq. ( [42] ) 
one can conclude that the p meson spectral function (^) used within the framework of the Borel QSR satisfies at the 
same time the following equation 

/°° ^T+i? = [n P (Q 2 ; o) + e(n Q1 (q 2 ; o) - n p (Q 2 ; o))] , (43) 

so that the in-medium modification of the p meson up to order T 2 stems from the admixture of the vacuum axial- 
vector correlator, which is not considerable as long as the temperature is sufficiently small in comparison with chiral 
transition temperature. 

The general trends of the p meson mass dependence on the density in the region of low temperature, which is relevant 
for the heavy- ion experiments at HADES, can be understood by means of the following approximation relying on the 
pole approximation of p? es 

F p e-^ 2 « Cb(l - e"^ 2 )M 2 - + ^ + ^ (44) 

where the four-quark condensate q^, the non-scalar condensate a 2 and the Landau damping term /f^ catt are given by 
eqs. ( Ek| [3l| , |39| ). Since the gluon condensate is weakly dependent on the density, and the quark condensate q2 and the 
non-scalar condensate 04 have rather small contributions, we have omitted the corresponding terms in eq. (|44[). It is 
clear from eq. ( f44| ) that in general the scattering term tends to increase the mass m p , while the non-scalar condensate 
02 tends to decrease the value of m p . When the value of the Borel mass M 2 is about in the center of the Borel 
window, i.e. M ~ 1 GeV, the numerical value of p^att IS ver y close to a2/M 2 , so that there is a cancellation of these 
terms in a wide range of nucleon densities. Therefore, the density dependence of m p is determined by the behavior 
of the four-quark condensate, in spite of the fact that the term qA c /{2M 2 ) is much smaller than p^ tt or a-ifM 2 
separately. Actually just the linear dependence of the four-quark condensate on the density leads to the approximate 
linear decease of m p with increasing density, as exhibited in fig. [l]b. 

Due to the crucial role of the four-quark condensate for the density dependence of m p it is useful to vary the 
parameters ko and kn which determine the numerical value of q^ and its dependence on the density (see eq. (^o|)). 
As seen in fig. ^| the main trend in decreasing the mass m p , obtained above for the parameterization kn/kq = 0.41 
||, is stable with respect to the variation of kq and kn within the reasonable numerical limits 2.0 < < 3 and 
< kat < 3. (Note that the value of kq is usually adjusted to the vacuum mass of the p meson and the ratio k^/ko is 
restricted by the conditions 54 < and 94 — > with increasing density and temperature.) In accordance with eq. ( |30| ) 
for <u and eq. (Q) the increase of the ratio kn/k-o leads to a stronger shift of m p at fixed density, as displayed in 
fig. [2| This further confirms the crucial role of the four-quark condensate for the in-medium behavior of the p meson 
mass. 

In a sufficiently dense nucleon medium one can expect a strong broadening of the p meson width due to inelastic 
pN scattering P,|l5|]. In principle, the inelastic processes should be included in the hadronic spectral function (p^), 
but in this case the OPE side should also include the higher order effects related, in particular, to two-pion states. 
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Since such effects are beyond the often used dilute gas approximation we do not concentrate here on the in-medium 
modifications of the vector meson widths and insist to postpone this to a separate and self-consistent study with 
higher order effects in the OPE side for finite density and temperature. 

The results of our QSR evaluation of the density and temperature dependence of the u> meson mass are exhibited 
in figs. ||a and b. Due to the comparatively large scattering term p^ tt 3> p^tt the density and temperature behavior 
of the lu meson mass differs essentially from the corresponding in- medium modifications of the p meson. In particular, 
for the parameterization of the four-quark condensate with kn/kq = 0.41, which is used above for the p meson, the 
mass of the u> meson increases with increasing density, as shown in fig. ||b for a wide region of the temperature. This 
can be still understood within the approximation (144), but now one has to take into account that the numerically 
large term p^f tt overwhelms the contributions from the non-scalar condensate a-2- A similar behavior of the ui meson 
mass was obtained in for the case T = 0. 

We consider also the sensitivity of the density dependence of the w meson with respect to various parameterizations 
of the g4 term. In fig. ^ we plot the density dependence of the uj meson mass at fixed temperature T — 20 MeV for 
various values of the parameter kn which reflects the strength of the density dependence of the four-quark condensate. 
As seen in fig. |] the in- medium modification of the u> meson changes even qualitatively under variation of the parameter 

kn (at fixed kq = 3.0): for kn ~ 2 the to meson mass increases with increasing density, while for kn > 2 it drops. 
Such a drastic change of the density dependence happens since the behavior of the density dependent part of the 
combination — p" c ^ tt + (I2/M 2 + q^f (2M 4 ) in the QSR eq. (Q) is governed by the parameter k-^. The obtained strong 
sensitivity of m u on the four-quark condensate holds also for higher temperature. From the above considerations 
one can conclude that the sign of the in-medium ui mass shift is directly related to the density dependence of the 
four-quark condensate. 

In contrast to the p and ui mesons the density dependence of the 4> meson is governed mainly by the in-medium 
chiral condensate of strange quarks m s (ss) llN ,Ti which has the dominant contribution in the OPE. As exhibited in 
figs. |^a and b, drops with increasing nucleon density. (To adjust the vacuum value of the (f> meson mass we employ 
kq = 0.6 and the values of parameters as listed in table 1.) This is in agreement with results of j|] obtained at T = 0. 
Due to the weak dependence of the chiral quark condensate on temperature the linear decrease of with increasing 
nucleon density holds in a wide range of temperatures up to 140 MeV. We also find a very tiny change of by 
variations of the four-quark condensate due to the density dependence. From the above analyses one can conclude 
that just the measurement of the <f> meson mass shift in a nucleon medium is most appropriate to search for in-medium 
modifications of the chiral quark condensate. In addition, the <f> meson is expected to keep its quasi-particle character 
as narrow resonance |l5|. Such measurements can be accomplished at HADES in heavy-ion experiments along with 
pion or proton induced (f> production at nuclei. 

VI. SUMMARY 

In summary we present for the first time a systematic analysis of the QCD Borel sum rules for light vector mesons at 
finite nucleon density and temperature. Our approach is based on the dilute gas approximation for a nucleon medium 
which is appropriate at not too large nucleon density (say, n < 2rt ) and temperature (say, T < 100 MeV). Such 
conditions are expected to be reached in heavy- ion experiments at HADES. We find that in- medium modifications 
of the p and uj mesons masses are dominated, within the Borel QSR approach, by the dependence of the four-quark 
condensate {{qq) 2 ) M „ t on density and temperature. In particular, the numerical value of the parameter k^v, which 
describes the strength of the linear density dependence of ((qq) 2 ) flnt T governs the decrease of the p meson as a function 
on the density. At the same time an experimental identification of the p meson mass shift is closely related to the 
problem of the strong broadening of its in-medium width. This needs separate considerations with the Borel QSR 
including higher density and temperature effects into the OPE side. (For a possible strategy to get access to a strongly 
broadened p meson by means of the double differential e + e~ spectra we refer the interested reader to [p4|.) 

For the u> meson the sign of the in- medium mass shift is changed by variation of the parameter k^. This also 
points to the crucial role of the poorly known four-quark condensate and makes definite predictions difficult. A direct 
measurement of the in-medium mass shift can give, within the considered framework, an important information on 
the in-medium behavior of the four-quark condensate. Since the difference of the vector and axial-vector correlation 
functions is proportional to the four-quark condensate the measured sign of the vector meson mass shift can serve, 
via the four-quark condensate, as a tool for determining how fast the nucleon system approaches the chiral symmetry 
restoration with increasing density. 

The in-medium <f> meson mass shift is directly related to the behavior of the chiral strange quark condensate 
^ns(ss)f_t N} T which decreases with density almost independently of the temperature in a wide region of density and 
temperature. This offers the chance to measure the density dependence of the chiral quark condensate via an in- 
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medium mass shift of the <j> meson in heavy-ion experiments, supposed the strangeness factor y is not too small J25|. 
Since at nuclear saturation density the mass shift is already noticeable, one can also complementary search for a <f> 
meson mass shift in hadron - nucleus reactions at the HADES detector under suitable kinematical conditions. 
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parameter 


numerical value 


reference 


a 3 


0.38 




26 




m q 


0.0055 GeV 




14 
™ 


1 


m s 


0.130 GeV 




1 


f« 


0.093 GeV 




14 


1 


M N 


0.938 GeV 




26 


f 


Mn 


0.770 GeV 




5|_ 


(Jjv 


0.045 GeV 




il 


y 


0.22 




1 

u 


lj 


m n 


0.138 GeV 






(uu) = (dd) 


(-0.245 GeV) 3 




cf. )5| 


(ss) 


0.7{uu)o 




27 
27 






(0.33 GeV) 4 






Au+d 
■™-2,N 


1.02 




5 




AU + d 

"-4.N 


0.12 




5 




AU + d 
2,7T 


0.97 




12 
12 




AU + d 


0.255 






A a 


0.1 




5| 




0.004 




i. 


/i 2,ir 


0.08 




1C 


] 


A a 


0.008 




19 


| 



TABLE I. Set of parameters used. 
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FIG. 2. The mass parameter m p as a function of kjv for /to = 3.0 and 2.0 (solid and dashed curves, respectively) at njv = 0.15 
fm~ 3 and T = 20 MeV. 
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